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Abstract 

We present two things in this discussion. First, we develop and prove the theory of 
the extension of the YFS Monte Carlo approach to higher order SU2L x Ui radiative 
corrections to the analogous higher order QCD radiative corrections. Contact is made 
with other pioneering soft gluon resummation theories in the literature. Second, semi- 
r^] analytical results and preliminary explicit Monte Carlo data are presented for the specific 
^1 example of the processes pp ^ ti + n{G) + X aX FNAL energies, where G is a soft gluon 
^ • and the respective event generator, ttpl.O, features realistic, event-by-event simulation of 
E> I multiple, soft, finite pt gluon effects in which the infrared singularities are canceled to 
^ ', all orders in a^. We comment briefly on the implications of our results on the CDF/DO 
^ \ observations and on their possible applications to RHIC physics and to LHC physics. 
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1 Introduction 



The problem of soft gluon resummation is well known [|l],0 and some of its many phe- 
nomenological applications are also: 

• FNAL tt production cross section higher order corrections ( the current situation [0, 
^ has the experimental cross section 6.21*1]^'^ to be compared with a theoretical 
prediction |^ of 5.1 ±0.5) and the attendant soft gluon uncertainty in the extracted 
value of mt = 0.1743 ± 0.0051TeV, where g ~ 2 - 3GeV of the the latter error 
could be due to soft gluon uncertainties. 

• RHIC hard scattering polarized pp scattering processes. 

• n{G) production in the hard nucleon-nucleon scattering processes which participate 
in the nucleus-nucleus collisions at RHIC, where G denotes a soft QCD gluon. 

• bb and J/"^ production by hard processes at FNAL. 

• heavy vector boson production at FNAL and RHIC, etc. 

For the LHC/TESLA/LC, the requirements on the corresponding theoretical precisions 
will be even more demanding and the QCD soft n{G) MC exponentiation which we discuss 
in the following will be an important part of the necessary theory - YFS [0] exponentiated 
0{a^)L corrections realized on an event-by-event basis. 

The results which we present also will ultimately allow us to investigate from a dif- 
ferent perspective some of the outstanding theoretical issues in perturbative QCD: 

• Treatment of the 2i(f)+n(G) phase space, where we have in mind our usual exact 
treatment of this phase space to be compared with the somewhat different approach 
of Catani and Seymour [§], for example - here, f denotes a generic fermion. 

• Approaches to resummation itself |I],P-pT|. 



No-go theorems, such as the result of Ref. |]T2[ , which requires that for calculations 



beyond 0{as), all initial state parton masses must vanish. 

Clearly, the results which we present herein have indeed a large arena for further devel- 
opment and application. For definiteness, we will use the process in Fig. 1 , Q{pi)Q{qi) 
i{p2)t{Pi) + Gi{ki) ■ ■ -Gnikn), as proto-typical process, where we have written the kine- 
matics as it is illustrated in the figure. This process, which dominates tt production 
at FNAL, contains all of the theoretical issues that we must face at the parton level to 
establish YFS QCD soft exponentiation by MC methods - extension to other related 
processes will be immediate. For reference, let us also note that, in what follows, we use 
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Figure 1: The process QQ t + t + n[G). The four-momenta are indicated in the 
standard manner: is the four-momentum of the incoming Q, q2 is the four-momentum 
of the outgoing t, etc., and Q = u, d, s, c, b, G. 



the GPS conventions of Jadach et al. [jT^ for spinors and the attendant photon and gluon 
polarization vectors that follow therefrom: 



V2 u_^{k)u^{(3y 



u^{k)ru^{C) 
V2 u^^{k)u^{0' 



(1) 



so that all phase information is strictly known in our amplitudes. This means that, 
although we shall use the older EEX realization of YFS MC exponentiation as defined 
in Ref. 11141, the realization of our results via the the newer CEEX realization of YFS 



exponentiation in Ref. ||T^ is also possible and is in progress 



We organize our presentation as follows. In the next section, we review the application 
of YFS MC methods in the EW-QED case. In Sect. 3, we prove that we can extend this 
application to the QCD theory. In Sect. 4 we illustrate the QCD extension by applying 
it to tt production at FNAL. Sect. 5 contains our summary remarks. 



2 Review of YFS Theory: An Abelian Gauge Theory 
Example 

In this section, for pedagogical reasons, we review the application of YFS exponentiation 
to the prototypical Abelian gauge theory example of e^e~ — > ff + n{'j) so that we set the 
stage for our proof of the extension of this theory to the prototypical non- Abelian gauge 



theory example of qq — > tt + n{G) in QCD in the next section. Specifically, in Refs. ||T6|-[T8 
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we have shown that the process e+(pi)e (gi) f{P2)f{<l2) + n{'y){ki, ■ ■ ■ ,kn), in the 
renormahzation group improved YFS theory ||16[, is represented by 



2a Re B+2a B 



^^^j f ^'^y ^iv(pi+<1i-P2-ci2-y.^k,)+D 



^n! y n fco J (27r)4- 

n=0 ^ j=l 3 ' 

R (h h \ ^^P2c/^g2 

Pn\'^l) • • • ; l^n) n n 

where the real infrared function B and the virtual infrared function B are given in Refs. [0, 



17|, p!9H22[| , and where we note the usual connections 

2a B = 



D 



m 
k^ 



-iy-k 



for the standard YFS infrared emission factor 

Pi 



S{k) 



a 
4^ 



QfQ 



k)) 



(2) 



'f\pi-k q,-k 



(3) 



if Q/ is the electric charge of / in units of the positron charge. For example, the YFS hard 

m 



photon residuals (3i in 



0, 1, 2, are given in Ref. 



for BHLUMI 4.04 so that this 
latter event generator calculates the YFS exponentiated exact 0{a) and LL 0[a^) cross 
section for Bhabha scattering using a corresponding Monte Carlo realization of . In the 
remainder of this paper, we will develop and apply the analogous theoretical paradigm to 
the prototypical QCD higher order radiative corrections problem for qq ^ tt + n{G). 



3 Extension to non-Abelian Gauge Theories: Proof 

Specifically, in Refs. [0], we have analyzed how in the special case of Born level color 



exchange one applies the YFS theory to QCD by extending the respective YFS IR sin- 
gularity analysis to QCD to all orders in a^. Here, we will take a somewhat different 
approach. We want to focus on the YFS theory as a general re-arrangement of renor- 
malized perturbation theory based on its IR behavior, just as the renormahzation group 
is a general property of renormalized perturbation theory based on its UV(ultra-violet) 
behavior. Thus, here we will keep our arguments entirely general from the outset, so that 
it will be immediate that our result applies to any renormalized perturbation theory in 
which the cross section under study is finite. 
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Let the amplitude for the emission of n real gluons in our proto-typical subprocess, 
Q? + Q? ^ ^^'^'^ + ''^{G), where a, a, 7, and 7 are color indices, be represented by 

^^(f" = 5^M(;J"^ (4) 

M^'^ is the contribution to }A^'^^ from Feynman diagrams with (. virtual loops. Sym- 
metrization yields 



Mi") 



where this last equation defines symmetric function of its arguments arguments 

ki, ...,k£. A will be our infrared gluon regulator mass for IR singularities; n-dimensional 



regularization of the 't Hooft-Veltman type is also possible as we shall see. 



We now define the virtual IR emission factor Sqcoik) for a gluon of 4-momentum /c, 
for the /c ^ regime of the respective 4-dimensional loop integration as in such that 

where we have now introduced the Born level color exchange condition asa^d^^^a 
for definiteness. (Henceforth, when we refer to — > gluons we are always referring for 
virtual gluons to the corresponding regime of the 4-dimensional loop integration in the 
computation of M^"\) 

In Ref. we have calculated SQcoik) using the running quark masses to regulate 



its collinear mass singularities, for example; n-dimensional regularization of the 't Hooft- 
Veltman type is also possible for these mass singularities and we will also illustrate this 
presently. 

We stress that Sqcoik) has a freedom in it corresponding to the fact that any function 
ASQCD{k) which has the property that lim^^o ^^^'S'qcd(^)Po"'' — added to it. 

Since the virtual gluons in pf''^ are all on equal footing by the symmetry of this 
function, if we look at gluon for example, we may write , for ki — s> (0, 0, 0, 0) = O while 
the remaining fcj are fixed away from O, the representation 

pf^ = SQcoih) * pf\{h: ■ ■ ■ , + Pl{ki, ■ ■ ■ , h-i] h) (7) 

where the residual amplitude (3j{ki, - ■ ■ , fc^-i; ki) will now be taken as defined by this last 
equation. It has two nice properties: 

• it is symmetric in its first i — 1 arguments 

• the IR singularities for gluon i that are contained in Sqcoike) are no longer con- 
tained in it. 
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We do not at this point discuss whether there are any further remaining IR singu- 
larities for gluon I in ■ ■ ■ , fc^-i; ki). In an Abehan gauge theory hke QED, as has 
been shown by Yennie, Frautschi and Suura in Ref. [0], there would not be any further 
such singularities; for a non-Abelian gauge theory like QCD, this point requires further 
discussion and we will come back to this point presently. 

We rather now stress that if we apply the representation (0) again we may write 

pf^ = SQCD{ke)SQCD{ke^i) * ^^"2(^1' " " " > ^^-2) 
+SQCD{ki)(3j_i{ki, • • • , ke_2; ^£-1) 
+'S'qcd(^£-i)/31-i(^i, • • • , ki_2; ki) 

+0i{ki, ■■■ , ki_2; ki_i,ki), (8) 

where this last equation serves to define the function 0f{k\, • • • , k^-2] ki-i, ki). It has two nice 
properties: 

• it is symmetric in its first £ — 2 arguments and in its last two arguments ki 

• the infrared singularities for gluons £ — I and £ that are contained in SqcDiki-i) and 
SgcDiki) are no longer contained in it. 

Continuing in this way, with repeated application of we get finally the rigorous, exact 
rearrangement of the contributions to p^^^ as 

e 

pP = SgcDiki) ■ ■ ■ SQCD{h)Pl + E n SQCD{k,)(i\{h) + ... 

i=l jj^i 

+(5',ik,,--- (9) 
where the virtual gluon residuals (3l{k[, - ■ ■ , A;-) have two nice properties: 

• they are symmetric functions of their arguments 

• they do not contain any of the IR singularities which are contained in the product 
SQcnik'i)--- SQcoiK). 

Henceforth, we denote /?• as the function Pi for reasons of pedagogy. We can not stress too 
much that (^) is an exact rearrangement of the contributions of the Feynman diagrams which 

(n) 

contribute to , it involves no approximations. Here also it is important to note that we 
do no enter into the question of the absolute convergence of these Feynman diagrams from the 
standpoint of constructive field theory. Yennie, Frautschi and Suura have already stressed 
that Feynman diagrammatic perturbation theory is non-rigorous from this standpoint. What we 
claim is that the relationship between the YFS expansion and the usual perturbative Feynman 
diagrammatic expansion is itself rigorous even though neither of the two expansions themselves 
is rigorous. 
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Introducing Q into @ yields a representation similar to that of YFS, and we will call it a 
"YFS representation" , 



oo 

= exp{asBQCD) \ (10) 

3=0 



where we have defined 



as{Q)BQCD = I -773, I ■ •. SQCDik) (11) 



(fc2-A2 + ie) 
and 

-?^=^ /n^'/^.^fci'---'^.)- (12) 

We say that ( |lO|) is similar to the respective result of Yennie, Prautschi and Suura in Ref. 0] and 
is not identical to it because we have not proved that the functions Pi{ki, ...,ki) are completely 
free of virtual IR singularities. What have shown is that they do not contain the IR singularities 
in the product Sqc niki) • • • SqcDiki) so that mj"^ does not contain the virtual IR divergences 
generated by this product when it is integrated over the respective 4j-dimensional j-virtual gluon 
phase space. In an Abelian gauge theory, there are no other possible virtual IR divergences; in 
the non-Abelian gauge theory that we treat here, such additional IR divergences are possible; 
but, the result (|l^) does have an improved IR divergence structure over (Q) in that all of the IR 
singularities associated with Sqcoik) are explicitly removed from the sum over the virtual IR 

improved loop contributions m^"^ to all orders in as{Q)- 

Turning now to the analogous rearrangement of the real IR singularities in the differential 
cross section associated with the A^^"\ we first note that we may write this cross section as 
follows according the standard methods 

2a,ReBQCD f rpl. " 



m=l ^ "1 ' ' 1=1 



P^^\Pl,P2,Ql,Q2,ku--- ,kn)^-^^, (13) 



where we have defined 



oo 



p'^''\Pi,P2,Qi,Q2,k,,... ,kn)= WY.^'ff (14) 

color, spin j=0 

in the incoming qq cms system and we have absorbed the remaining kinematical factors 
for the initial state flux, spin and color averages into the normalization of the amplitudes 
_A/iW for reasons of pedagogy so that the p^^^ are averaged over initial spins and colors 
and summed over final spins and colors. We now proceed in complete analogy with the 
discussion of p^"^ above. 
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Specifically, for the functions p^''^\pi,p2, qi, q2, ki, - ■ ■ , kn) = p^"'\ki, ■ ■ ■ , kn) which are 
symmetric functions of their arguments /ci, ■ ■ ■ , kn-, we define first, for n = 1, 

hm k'(-p^'\k)-~SQCD{k)p^'^)=Q. (15) 

where the real infrared function SQcnik) is rigorously defined by this last equation and 
is explicitly computed in DeLaney et al. [Q; like its virtual counterpart SqcDik) it has a 
freedom in it in that any function ASqcnik) with the property that lim|g|_,Q k'^ASqcnik) = 
may be added to it without affecting the defining relation ([T5|) . 



As we show in Ref. [|^, we can again repeat the analogous arguments of Ref. to 
get the "YFS-hke" result 

dlTexp = ^ da"" 

n 

= eSUMiR(QCD) ^ ij^'^h. I J^^f,iyiPi+P2-Qi-Q2-T,kj)+DQCD ^^g^ 



* Pn{kl, • • • , kn) Q Q 
1^2 H 



2 



with 



SUMir{QCD) = 2a,ReBQCD + 2a.5QCD(i^max), 

/d^k ~ 
-I^SQCD{k)6{Kjna.^ - k), 

— 5QCD(fe) [e-*^-^ - e{K^,^ - k)\ , (17) 

i;3l = da^l-SQCD(^)d<7B, ... (18) 

where the are the QCD hard gluon residuals defined above; they are the non-Abelian analogs 
of the hard photon residuals defined by YFS. Here, for illustration, we have recorded the relation- 
ship between the n = 0, 1 through 0{as) and the exact one-loop and single bremsstrahlung 
cross sections, da^^~^°°^\ da^^ , respectively, where the latter may be taken from Nason et al. [26 
and Beenakker et al. We stress two things about the right-hand side of (Il6|) : 



It does not depend on the dummy parameter K^ax which has been introduced for can- 
cellation of the infrared divergences in SUMir{QCD) to all orders in as{Q) where Q is 
the hard scale in the parton scattering process under study here. 
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Its analog can aslo be derived in our new CEEX \14\ format 



We now return to the property of (16) that distinguishes it from the Abelian result derived 
by YFS - namely, the fact that, owing to its non- Abelian gauge theory origins, it may very well 
be that there are infrared divergences in the which were not removed into the Sqcd, Sqcd 
when these infrared functions were isolated in our derivation of ([TBI). 

More precisely, the left-hand side of (|TBp is the fundamental reduced parton cross 
section and it should be infrared finite or else the entire QCD parton model has to be 
abandoned. 

There is an observation in the literature [1^ that unless we use some collinear mass 
singularity regulator other than the incoming light quark running masses, the reduced par- 
ton cross section on the left-hand side of (|16]) diverges in the infrared regime at 0{al{Q)). 
We do not go into this issue here but use n-dimensional methods to regulate such diver- 
gences while setting the quark masses to zero as that is an excellent approximation for 
the light quarks at FNAL energies - we take this issue up elsewhere. 

From the infrared finiteness of the left-hand side of (|16]) and the infrared finiteness of 
SUMjji{QCD), it follows that the quantity 

c/o-exp = exp[SUMiR(QCD)]c?(Texp 

must also be infrared finite to all orders in a,. 



As we assume the QCD theory makes sense in some neighborhood of the origin for 
as, we conclude that each order in must make an infrared finite contribution to da cxp- 
At , the only contribution to da^xp is the respective Born cross section given 

by in (|16D and it is obviously infrared finite, where we use henceforth the notation 
Pn^ to denote the 0{al{Q)) part of /3„. Thus, we conclude that the lowest hard gluon 
residual is infrared finite. 

Let us now define the left-over non- Abelian infrared divergence part of each contribu- 
tion /3n^ via 

where the new function (3^ is now completely free of any infrared divergences and the 
function Dj3n'^ contains all left-over infrared divergences in l3n' which are of non-Abelian 
origin and is normalized to vanish in the Abelian limit fabc ~^ where fahc are the group 
structure constants. 

Further, we define Dfin'^ by a minimal subtraction of the respective IR divergences 
in it so that it only contains the actual pole and transcendental constants, 1/e — Ce for 
e = 2 — d/2, where d is the dimension of space-time, in dimensional regularization or In 
in the gluon mass regularization. Here, Ce is Euler's constant. 
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For definiteness, we write this out explicitly so that there can be no confusion about 
what we mean 



/n+e 
i=i 



where the coefficient functions c/"'^ are independent of A for A ^ and dPh is the 
respective n-gluon Lorentz invariant phase space. 

At 0(q;"((5)), the IR finiteness of the contribution to rfaexp then requires the contri- 
bution 



i=0 j = t^j^J^max 3 j = 3=1+1 



Pr ,ke+,)^-^ (19) 



to be finite. 

From this it follows that 

„ n £ „ ,3, n—£ ^ i+i 

1=0 j = l-^ f'j^J^max J j=0 j=£+l 

-TT^DPl '{ke+i,... ,ke+i)—^—^ (20) 



is finite. Since the integration region for the final particles is arbitrary, the independent 
powers of the IR regulator In(A^) in this last equation must give vanishing contributions. 

This means that we can drop the D/5n^ from our result (|I6|) because they do not make a 
net contribution to the final parton cross section (Jexp- We thus finally arrive at the new 
rigorous result 

daexp = ^ da"- 



00 



n=0'' j=l ^^^-^ 

~- ,.d^P2d^Q2 
* P^[Ki, . . . ,Kn) Q Q 

^2 

where now the hard gluon residuals P^iki, . . . , kn) defined by 

Pni^lj ■ ■ ■ 5 — / ^ Pn (^1' • • • i ^n) 
£=0 

are free of all infrared divergences to all orders in as{Q). This is a basic result of this paper. 
We note here that, contrary to what we claimed in the Appendix of the first paper in Refs. [p4|. 
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the arguments in Refs. [24| are not sufficient to derive the respective analog of eq.(^); for, they 
did not reahy expose the compensation between the left over genuine nonAbelian IR virtual and 
real singularities between J dPh(3n and / dPhf3n+i respectively that really distinguishes QCD 
from QED, where no such compensation occurs in the /3„ residuals for QED. Further discussion 
of the many implications of eq. (pll) is given elsewhere 



4 YFS Exponentiated QCD Corrections to tt Pro- 
duction at High Energies 

We shall realize the result above as it is applied to the process in Fig. 1 at high energies by the 
Monte Carlo event generator methods of Ref. p^,p^. The first application of these methods 
to QCD processes has already been reported in DeLaney et al. ||2^. Here, we shall apply the 
general result in the latter reference to the tt production at the parton level and show how to 
synthesize that parton level result with the realistic parton distributions to arrive at an event 
generator for pp ^ ti + X. Similar results will hold for pp incoming states. Sample MC data 
will be illustrated. We refer to the respective MC event generator as ttpl.O. It is in the EEX 
format but a CEEX version is imminent. It will be available from the authors soon. 

The starting point in our MC realization of basic result for tt production at the hadron level 



is its realization for the respective parton level processes. In Ref. [24| we have shown how to 
construct such parton level event generators for the processes such as qq ti+n{G) where q is 
a light quark. Here, we extend these results to the process G + G — > ti+n{G). This is presented 
in detail elsewhere mM- 



Following the procedures in Ref. [24|, we then use the Monte Carlo algorithm presented 
in S. Jadach et al. |17, 1^ to realize the result derived in the previous section for both the 
qq ^ tt + n{G) and GO ti + n{G) subprocesses on an event-by-event basis in which infrared 
singularities are now canceled to all orders in q^. 

In order to apply these parton level results to the desired hadron level cross section a{pp — > 
ti + X), we use the standard formula 

a{pp ^ti+X) = Fi{xi)Fj{xj)da'^^^ ijdxidxj (22) 

where Fi{Fj) is the structure function of parton in p(p) and where i j result 
derived above for the tt production subprocess with the incoming parton-i,parton-j initial state 



when the DGLAP synthesization procedure presented in Ref. |28] is applied to it to avoid over- 
counting resummation effects already included in the structure function DGLAP evolution. In 
operational terms, the DGLAP synthesized result o^^^p j^- is obtained from the result o'cxp/y 

of 

the previous section by substituting the functions Sqqj^, Bq^,j^ for Sqcd, Bqcd, respectively, 
where the former functions are given in Ref. and do not contain big logs already contained 
in the structure functions, for example. 



The formula in ( p2D we have realized by MC methods by extending the MC realizations of 
^ ) for the qq and GG subprocesses to the respective realizations of da'^^^ , ij = qq, GG in 
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the presence of the additional two-dimensional structure function distribution in a standard way 
which is already illustrated for example in Ref. |^]. In this way, we arrive at the MC event 
generator ttpl.O, which features YFS-style exponentiated multiple soft gluon radiative effects in 
pp— >tt-|-X on an event-by-event basis in which infrared singularities are canceled to all orders 
in as- We now illustrate its application with both semi-analytical results and with with some 
preliminary sample MC data at FNAL energies. 

As a first check of the approach, we have realized the formula in (^Tj) by semi-analytical 
methods. This is useful for a number of reasons. For example, it allows us to check the normal- 
ization of our work with that of other resummed calculations in the literature. An important 
outcome is an estimate of the size of the 0(Q^),n > 2 contribution to a{tt) at FNAL B. 



Specifically, as we show in Ref. |3C|, our semi-analytical result for the ratio of the YFS 
exponentiated cross section to the corresponding Born cross section is given by 



nls 
exp 



exp 



[{2Cf 



IT 

1.086, as 

1.103, as 

[1.110, as 



asiVs), 

as{2mt), 

as{mt). 



1 



Cf] 



(23) 



This implies that the prediction for the contribution to (T{pp — > tt) from {0{a^),n > 2) is 



,nls 
exp 



a 



vr 



1-^[(2C,,--Ca)- 



vr 



- Cij ])daB{xiXjs)dxidx 



J' 



(24) 



where 



Cf, ij = qq, qq, 



Ca, 



GG. 



(25) 



From this we get that (0(a"),n > 2) contributes .006-. 008 of the NLO cross section, in 
agreement with Catani et al. in Ref. 0. 

We have also generated preliminary MC data using the proto-typical version of ttpl.O. 
A detail presentation of such data will appear elsewhere [|T5|. Here, in view of its impor- 
tance, we show preliminary data for the px spectrum from ttpl.O. This is shown in Fig. 
2. 
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Figure 2. n{G) transverse momentum distribution for in pp^ tt +X+n{G). The result 
is a preliminary. 

What we see is that, indeed, the soft n{G) effects do modulate significantly the pt spec- 
trum on the scale of a few GeV - they must be taken into account properly to understand 
observables derived from them on this level. 

5 Conclusions 

We conclude that YFS theory, both the EEX and CEEX formulations, extends to QCD. In 
our work, a full MC event generator realization of this extension is near. Our early semi- 
analytical results agree with the literature on tt production at FNAL energies. Preliminary 
MC data show that the n{G) pt is significant in the latter process. We are currently 
pursuing the attendant ramifications for LHC and RHIC as well. 
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Figure Captions 



Figure 1. The process QQ t + t + n{G). The four-momenta are indicated in the 
standard manner: qi is the four-momentum of the incoming Q, q2 is the four-momentum 
of the outgoing t, etc., and Q = u,d, s, c, b, G. 

Figure 2. n{G) transverse momentum distribution for in pp^ tt +X+n(G). The result 
is a prehminary. 
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